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NOTATION 


Van  Driest' s  damping  factor,  A  =  26v 

Length  of  major  elliptical  axis  at  given  x/L 

^  2 

Turbulence  structure  parameter,  a^  =  -u'v'/q 

Effective  displacement  thickness,  see  Equation  (3) 

Length  of  minor  elliptical  axis  at  given  x/L 

Pressure  coefficient,  =  (p-pQ)/(l/2pU^)  =  1  -  (Ue/UQ)^ 

Metric  coefficients 

Geodesic  curvatures  of  the  curves  z  =  constant  and 
x  =  constant,  respectively 

Functions  of  the  geodesic  curvatures  and  metric  coefficients 
Total  body  length 

Mixing  length  parameter:  In  the  inner  region 
I  *  0,4  y  [ 1-exp (-y /A) } 

In  the  outer  region 


Coordinate  measured  normal  to  the  body  profile  in  the  y-z  plane 

Measured  local  static  pressure 
Measured  ambient  pressure 

Measured  static  pressure 

Measured  dynamic  total  pressure 
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Turbulence  parameter,  q  =  u^  +  v'  +  v' 
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Eddy  viscosity  in  the  inner  and  outer  regions,  respectively. 
Equation  (2) 

Angular  coordinate  measure  in  the  y-z  plane  from  the  z— axis 
line  joining  the  surface  offset  and  elliptic  center 

Angle  between  the  x  and  z  coordinates 

Effective  displacement  area 

Kinematic  viscosity  of  the  fluid 

Mass  density  of  the  fluid 
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ABSTRACT 


A  comprehensive  set  of  experimental  pressure,  velocity,  and 
turbulence  data  are  presented  across  the  stern  of  a  three- 
dimensional  model  having  3:1  elliptic  transverse  cross  sections. 

The  axisymmetric  displacement  body  concept  is  extended  to  three- 
dimensions  and  the  pressure  and  velocity  data  are  compared  with 
the  predictions  of  existing  three-dimensional  theoretical  methods. 
The  surface  pressures  for  the  displacement  body  are  found  to  model, 
satisfactorily,  the  measured  pressure  coefficients  in  all  regions 
except  over  the  aft  7  percent  of  body  length.  In  this  tail  region, 
the  boundary  layer  is  much  thicker  than  the  cross  section  dimensions 
and  the  theory  overpredicts  the  measured  distributions  of  the  mean 
velocity.  Agreement  is  particularly  poor  in  the  inner  region  of 
the  tail  boundary  layer,  indicating  a  need  to  examine  the  eddy 
viscosity  model  currently  used  in  computing  the  thick  stern  boundary 
layer  of  three-dimensional  models.  As  was  found  in  the  axisymmetric 
case,  the  measured  values  of  turbulence  intensity,  eddy  viscosity, 
and  mixing-length  parameters  in  the  stern  region  are  much  smaller 
than  those  of  a  thin  boundary  layer. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  funded  under  the  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center's  Independent  Research  Program,  Program  Element 
61152N,  Project  Number  ZR  000  01,  and  Work  Unit  1542-103. 

INTRODUCTION 

Many  single-screw  ship  propellers  operate  inside  of  thick  stern  boundary  layers. 

Satisfactory  predictions  of  turbulent  boundary-layer  characteristics  can  be  made  for 

the  forward  portions  of  a  body  by  solving  the  boundary-layer  equations  in  either 

integral  or  differential  forms.  However,  at  the  ship  stern,  the  thickness  of  the 

boundary  layer  increases  rapidly,  mainly  due  to  the  diminishing  cross-sectional  area. 

The  thickness  of  the  stern  boundary  layer  usually  exceeds  the  thickness  of  the  body. 

Detailed  measurements  of  the  turbulent  boundary-layer  characteristics  in  the  thick 

*1  2 

stern  boundary  layers  of  axisymmetric  bodies  have  been  made  by  Huang  et  al.  ’  in 

order  to  gain  insight  into  the  physics  of  thick  stern  boundary  layers.  These 

measurements  have  been  used  to  validate  the  displacement  body  concept  as  suggested 
3  4 

by  Preston  and  Lighthill  for  solving  viscid-inviscid  flow  interaction  and  an 
improved  turbulence  model  has  been  obtained  for  computing  thick  axisymmetric  boundary 


*A  complete  listing  of  references  is  given  on  page  105. 
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layers  on  two  convex  sterns  and  one  concave  stern.  ’  The  present  work  is  an  initial 

investigation  into  extending  to  three-dimensions  the  previous  studies  on  axisymmetr ic 

1  2 

bodies  by  Huang  et  al.  ’ 

Experiments  have  been  made  to  measure  the  flow  across  the  thick  stern  boundary 
layer  of  a  three-dimensional  body  having  a  3:1  elliptical  transverse  cross  section. 

A  10.06  ft  (3.07  m)  fiberglass  model  was  tested  in  the  Center's  Anechoic  Flow 
Facility  at  a  speed  of  100  ft/sec  (30.48  m/s),  resulting  in  an  overall  Reynolds 
number  based  on  length  of  6.5  x  10^.  Pressure  taps,  embedded  in  the  model,  were 
used  to  measure  the  pressure  distribution  on  the  surface.  Velocity  and  turbulence 
characteristics  were  measured  using  a  two-element  hot-film  sensor  and  were  analyzed 
with  an  on-line  computer.  Measurements  include  mean  velocity  profiles,  turbulence 
intensities,  Reynolds  stresses,  eddy  viscosity,  and  mixing  length. 

Several  experimental  quantities  are  compared  with  data  from  existing  theoreti¬ 
cal  methods  using  an  iterative  scheme.  The  potential  flow  distribution  on  the  body 
surface  is  computed  using  the  XYZ  Potential  Flow  (XYZPF)  computer  code  of  Dawson  and 

Dean."*  An  initial  boundary-layer  computation,  using  the  McDonnell  Douglas 
6  2 

Corporation,  Cebeci,  Chang,  Kaups  (C  K)  computer  code,  is  made  using  the  potential- 
flow  pressure  distribution  on  the  body.  Flow  separation  is  predicted  for  this  model 

o 

by  the  C  k  code  at  axial  locations  greater  than  4  percent  of  the  body  length  and 
angular  locations  greater  than  75  degrees.  Excessive  boundary-layer  growth  in  the 
separated  region  caused  the  boundary-layer  calculation  to  abort  prematurely  at  81 
percent  of  the  body  length.  Predictions  of  the  effective  displacement  thickness  for 
the  remaining  19  percent  of  the  body  length  are  obtained  by  extrapolation.  The 
potential  and  boundary-layer  flow  calculations  are  repeated  once  for  a  modified  body 
and  wake  geometry,  formed  by  adding  the  computed  effective  displacement  thickness. 
Comparison  of  predicted  and  measured  results  shows  that  this  procedure  predicts 
accurate  values  of  pressure  over  the  forward  93  percent  of  the  body  and  accurate  mean 
velocity  profiles  in  locations  where  the  boundary  layer  is  thin  compared  with  cross- 
sectional  area.  The  measured  eddy  viscosity  distribution  is  compared  with  the  thin 
boundary— layer  model  of  Cebeci  '  and  is  found  to  be  smaller  than  predictions. 

In  the  following  sections,  the  experimental  techniques  and  model  geometry  are 
given  in  detail.  The  experimental  data  are  presented  and  compared  with  theoretical 
predictions.  The  raw  data  and  derived  results  are  given  in  tabular  form  for  inde¬ 
pendent  use  by  other  investigators. 
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WIND  TUNNEL  AND  MODEL 


The  experimental  investigation  was  conducted  in  the  DTNSRDC  Anechoic  Wind  Tunnel 
Facility.  The  wind  tunnel  has  a  closed  jet  test  section  that  is  8  ft  (2. A  m)  square 
and  13.75  ft  (A.  19  m)  long.  The  corners  have  fillets  which  are  carried  through  the 
contraction.  The  test  section  is  followed  by  an  acoustically-lined  large  chamber 

23.5  ft  (7.16  m)  long.  It  was  found  previously,  by  Huang  et  al.,^  that  the  ambient 

free-stream  turbulence  levels,  ^u'2/L’q jx  100,  are  0.075,  0.090,  0.100  and  from 
0.12  to  0.15  for  free-stream  velocities,  U  ,  of  2A.A,  30.5,  38.1,  and  A5.7  m/s, 
respectively.  Integration  of  the  measured  noise  spectrum  levels  in  the  test  section 

from  10  to  10,000  Hz  indicated  that  the  typical  background  acoustic  noise  levels  at 

2 

30.5  m/s  were  about  93  dB  re  0.0002  dyne/cm  (0.0002  Pa).  These  levels  of  ambient 
turbulence  and  acoustic  noise  were  considered  low  enough  so  as  not  to  unfavorably 
affect  the  measurement  of  boundary-layer  characteristics.  The  maximum  air  speed  that 
can  be  achieved  is  200  ft/sec  (61  m/s);  in  the  present  experiments  the  wind  tunnel 
velocity  was  held  constant  at  100  ft/sec  (30. A8  m/s). 

A  simple  three-dimensional  body,  having  a  3:1  elliptic  transverse  cross  section 
with  a  bow  entrance  length  of  6.23  ft  (1.897  m) ,  was  used  for  the  present  experi¬ 
mental  investigation.  The  total  model  length  is  10.06  ft  (3.07  m)  with  a  maximum 
major  axis  of  1.588  ft  (0.A8  m)  and  a  maximum  minor  axis  of  6.35  in.  (16.12  cm). 

A  schematic  of  the  three-dimensional  afterbody  with  the  3:1  elliptic  cross  section 
is  shown  in  Figure  1.  The  major  and  minor  elliptic  axes  are  shown  in  Figure  1  as 
a  and  b,  respectively.  The  model  is  shown  in  the  anechoic  wind  tunnel  facility  in 
Figure  2.  The  support  struts  shown  in  the  figure  are  not  the  struts  used  for  this 
experiment.  Model  offsets  are  presented  in  Table  1. 

The  model  was  supported  by  two  streamlined  struts  separated  by  one-third  of  the 
model  length.  The  struts  are  0.5-in.  (1.27-cm)  thick  with  a  1.5-in.  (3.81-cm)  chord 
upstream  and  2.25-in.  (5.72-cm)  thick  with  a  6.0-in.  (15.2A-cm)  chord  downstream. 

The  model  is  designed  to  rotate  90  degrees  radially  about  a  center  axis  to  permit 
vertical  traversing  normal  to  the  surface  pressure  taps  (see  section  on 
Instrumentation).  The  disturbances  generated  by  the  supporting  struts  were  within 
the  region  below  the  horizontal  centerplane.  Therefore,  all  of  the  experimental 
data  were  taken  above  the  model  on  the  vertical  centerplane  along  the  upper  meridian 
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where  there  was  no  effect  from  the  supporting  struts.  One-half  of  the  model  length 
protruded  beyond  the  closed  jet  working  section  into  the  open-jet  section.  The 
ambient  static  pressure  coefficients  across  and  along  the  entire  open-jet  chamber 
(7.2  m  x  7.2  m  x  6.4  m)  were  found  to  vary  less  than  0.3  percent  of  the  dynamic 
pressure.  Tunnel  blockage  and  longitudinal  pressure  gradient  effects  along  the 
tunnel  length  were  almost  completely  removed  by  testing  the  afterbody  in  the  open- 
jet  section. 

The  location  of  the  boundary- layer  transition  from  laminar  to  turbulent  flow  was 

artificially  induced  by  a  0.024-in.  (0.61-mm)  diameter  trip  wire  located  at  x/L  = 

0.05.  Huang  et  al.  found  that  the  trip  wire  effectively  moved  the  location  of  the 

virtual  origin  to  x/L  =  0.015  for  axisymmetric  models  at  a  length  Reynolds  number 
6  8 

of  5.9  x  io  .  The  virtual  origin  for  the  turbulent  flow  is  defined  such  that  the 

sum  of  the  laminar  frictional  drag  from  the  nose  to  the  trip  wire,  the  parasitic 

drag  of  the  trip  wire,  and  the  turbulent  frictional  drag  aft  of  the  trip  wire  is 

equal  to  the  sum  of  the  laminar  frictional  drag  from  the  nose  to  the  virtual  origin 

and  the  turbulent  frictional  drag  from  the  virtual  origin  to  the  after  end  of  the 

model.  The  virtual  origin  locations  for  the  three-dimensional  body  are  expected  to 

be  different  for  different  streamlines.  Due  to  the  limited  number  of  grid  locations 

2 

used  in  the  present  calculation,  the  location  of  the  transition  for  the  C  K  boundary- 

layer  calculation  is  set  at  a  constant  value  of  x/L  =  0.030.  The  computed  differ- 

1  2 

ences  in  velocities  using  x/L  =  0.01  and  x/L  =  0.03,  for  axisymmetric  body  1, 
are  found  to  be  less  than  0. 1  percent  of  the  free-stream  velocities  in  the  tail 
region.  Thus,  the  error  of  using  the  constant  transition  location  of  x/L  =  0.03 
for  the  present  C  K  computation  is  expected  to  be  negligible. 

INSTRUMENTATION 

A  series  of  0.031-in.  (0.8-mm)  diameter  pressure  taps  were  embedded  normal  to 
the  surface  of  the  stern  at  nine  x/L  locations.  When  the  model  was  rotated  about 
its  axis,  the  pressure  taps  were  at  the  upper  meridian  location.  Additional  taps 
were  added  for  model  alinement;  see  Figures  3  and  4.  The  model  was  alined  by 
balancing  the  surface  static  pressure  about  a  line  of  symmetry.  From  Figure  3,  the 
model  is  alined  when  symmetrically  located  pressure  taps  at  c  and  d,  and  at  e  and  f, 
give  equal  pressures,  i.e.,  p  (c)  =  p  (d),  p  (e)  =  p  (f).  The  model  was  rotated  to 
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eight  test  positions  and  the  alinement  was  checked  by  the  pressure  balance  technique. 

A  Preston  tube  using  a  0.072-in.  (1.83-mm)  inside  diameter  was  attached  and  alined 

with  the  flow  at  the  pressure  taps  to  measure  the  shear  stress.  The  Preston  tube 

was  calibrated  in  a  1-in.  (2.54-cm)  diameter  water-pipe  flow  facility  described 

9 

by  Huang  and  von  Kerczek.  These  pressure  taps  were  connected  to  a  multiple  pres¬ 
sure  scanivalve  system  that  takes  one  integral  pressure  transducer  with  its  zeroing 
circuit  and  measures  a  single  pressure  in  sequence  along  the  stern  upper  meridian. 

The  pressure  transducer  was  designed  for  measuring  low  pressures  of  up  to  1  psi 

_3 

(6.895  x  10  Pa).  The  zero-drift  linearity,  scanivalve  hysteresis,  and  pressure 
transducer  zeroing  circuit  were  carefully  checked  and  the  overall  accuracy  was  found 
to  be  within  0.5  percent  of  the  dynamic  pressure. 

The  mean  axial  and  radial  velocities  and  the  turbulence  intensities  for  the 
Reynolds  stress  calculations  were  measured  by  a  TSI,  Inc.  Model  1241-20  "X"  type 
hot-film  probe.  The  probe  elements  are  0.002  in.  (0.05  mm)  in  diameter  with  a 
sensing  length  of  0.04  in.  (1.0  mm).  The  spacing  between  the  two  cross  elements  is 
0.04  in.  (1.0  mm).  A  typical  schematic  of  the  hot-film  probe  used  is  shown  in 
Figure  5.  A  two-channel  hot-wire  and  hot-film  anemometer  with  linearizers  was  used 
to  monitor  the  response  of  the  hot-film  probe.  A  temperature  compensating  sensor 
(probe)  was  used  with  each  hot-film  element  to  regulate  the  operating  temperature  of 
the  sensor  with  changes  in  air  temperature.  The  "X"  hot  film  and  its  temperature- 
compensated  sensor  were  calibrated  together  through  the  expected  air  temperature- 
range  and  supplied  with  their  individual  linearization  polynomial  coefficients  at 
the  factory. 

The  frequency  response  of  the  anemometer  system,  for  reliable  measurements 
claimed  by  the  manufacturer,  is  0  to  100  kHz.  Calibration  of  the  "X"  hot  film  was 
made  before  and  after  each  set  of  measurements.  It  was  found  that  the  hot-film 
anemometer  system  had  a  +0.5  percent  accuracy,  +0.75  ft/sec  (+0.23  m/s)  accuracy 
at  the  free-stream  velocity  of  150  ft/sec  (45.72  m/s),  during  the  entire  experiment. 
An  estimate  was  made  of  the  crossflow  velocity  by  yawing  the  "X"  hot-film  probe  in 
the  free  stream.  It  was  found  that  the  crossflow  velocities  were  about  one  percent 
of  the  free-stream  velocity. 

The  linearized  signals  were  fed  into  a  Time/Data  Model  1923-C  real-time 
analyzer.  Both  channels  of  the  analog  signal  were  digitized  at  a  rate  of  128  points 
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per  second  for  8  sec.  These  data  were  immediately  analyzed  by  a  computer  to  obtain 
the  individual  components  of  mean  velocity,  turbulence  fluctuation,  and  Reynolds 
stress  on  a  real  time  basis. 


A  traversing  system  with  a  streamlined  strut  was  mounted  on  a  guide  plate  that 
permitted  the  traverse  to  be  locked  in  various  stationary  positions  parallel  to  the 
longitudinal  model  axis. 


DISPLACEMENT  BODY  METHOD 

The  theoretical  method  evaluated  in  this  report  is  an  initial  attempt  at  ex¬ 
tending  to  three-dimensions  the  displacement  body  concept  described  by  Wang  and 
Huang*9  and  by  Huang  et  al.,*’^  for  axisymmetric  bodies.  The  pressure  distribution 
is  calculated  using  the  XYZ  Potential  Flow  (XYZPF)  computer  code  of  Dawson  and 
Dean.  The  input  offsets  to  the  XYZPF  code  are  given  in  Table  1.  The  boundary- 

layer  flow  over  the  body  is  calculated  by  using  the  differential  method  of  Cebeci, 

2  6 

Chang,  and  Kaups  (denoted  C  K) .  The  flow  in  the  wake  is  modeled  only  in  the  near 

wake  region  of  0.93  £  x/L  £  1.05. 

2 

The  C  K  method  consists  of  using  Keller's  two-point  finite  difference  method 
and  Cebeci  and  Stewartson's  procedure9  for  computing  flows  in  which  the  transverse 
velocity  component  contains  regions  of  reverse  flow  to  solve  three-dimensional 
boundary- layer  equations.  The  governing  equations  for  three-dimensional  incom¬ 
pressible  laminar  and  turbulent  flows  are  given  by 
Continuity  Equation 


h  (uhz sin  9)  +  h  <whi 


sin  9)  +  (vh^l^  sin  0)  =  0 


x-Momentum  Equation 


u  9u  ,  w_  9u  9u 
hj  9x  9z  V  9y 


cot  9  +  ^w 
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esc  9  + 


(la) 


csc^  9  8  (p/p)  cot  9  esc  9  9  (p/p)  9  /  9u  ■? 

hj  9x  h2  9z  9y  V  9y  u  V 


(lb) 


z-Momentum  Equation 


u  3w  ,  w_  3w  3w 
hj  3x  h£  3z  V  3y 


cot  0  +  K^u 
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esc  0  +  ^2^uw 


cot  0  esc  0  3 (p/p) 
h^  3x 


csc^  3, (p/p.)  +  a_  /  3w  _ 

h2  3z  3y  \  3y  / 


(lc) 


where  u,  v,  and  w  = 
x,  y,  and  z  = 


P 

P 


K12’  K21 


velocity  components  in  the  x,  y,  and  z  directions,  respectively 

nonorthogonal  boundary-layer  coordinates,  as  given  in 
Reference  6 

fluid  density 

pressure  on  the  body 

metric  coefficients 

geodesic  curvatures  of  the  curves  z  =  constant  and 
x  =  constant,  respectively 

functions  of  the  geodesic  curvatures  and  metric  coefficients 


0  =  angle  between  the  coordinates  x  and  z 

V  =  kinematic  viscosity  of  the  fluid 


-* - X  - X — X 

U  V  ,  V  w 


Reynolds  stresses 


The  eddy-viscosity  concept  is  used  to  relate  the  Reynolds  stresses  to  the  mean 
velocity  profiles  by 


- ? - T 

U  V 


£  3u 

°  3y 


inner  region  0  <_  y  £  y£ 


outer  region  y^  <  y 


(2) 


where 


1/2 


which  is  the  eddy  viscosity  in  the  inner  region 


EQ  *  0.0168  J  (ute~ut)dy 
0 

<  2.  2 q\  1/2 

a  =  (a  +w  +2u  w  cos  0) 
te  e  e  e  e 

=  (u^+w^+2uw  cos  0)*^ 

y  is  the  value  of  y  at  which  e.  =  e 
c  3  i  o 

The  displacement  model  presently  used  for  this  body  adds  the  theoretical  effec¬ 
tive  displacement  thickness  (defined  below)  to  the  body  surface  along  the  major  (y)- 
and  minor  (z)-axes  of  the  elliptic  cross  section.  The  surface  profile  along  each 
of  these  axes  is  extended  by  hand-fairing  from  the  location  of  separation,  or  93 
percent  of  body  length  (whichever  occurs  first),  to  5-percent  aft  of  the  body, 
resulting  in  an  open  body.  An  elliptical  cross  section  is  defined  between  the  off¬ 
sets  of  the  major  and  minor  axes. 

2 

The  C  K  computer  program  has  been  modified  to  compute  the  effective  displace¬ 
ment  thickness  a*  at  the  major  and  minor  axes  along  the  axial  length  of  the  body. 

The  definition  for  a*,  which  is  similar  to  the  axisymmetric  expression,  is 
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where  r  “  radius  of  curvature  at  the  particular  axis  of  interest  in  the  y~z  plane, 
C  i.e.. 


[»  (*)*] 


“y-axis 


[>♦  (U)2] 


'z-axis 


f°  /  ut  \ 

A*  =  effective  displacement  area.  A*  =J  ^1-  — — j rdy 
a  =  angle  between  the  body  surface  and  the  body  axis, 

a  3  tan_1  (s£)  or  01  =  tan_1  (sf ) 


r  =  r  +  y  cos  cx 
c 

y  =  normal  distance  from  the  wall 

Unlike  the  procedure  for  an  axi symmetric  body,  which  uses  an  iterative  procedure 
consisting  of  the  calculation  of  pressure  and  boundary-layer  flow  over  successive 
displacement  bodies,  the  present  scheme  for  three-dimensional  bodies  uses  only  one 
iteration.  The  uncertainties  in  defining  the  displacement  body  in  the  region  between 
the  major  and  minor  axes  and  in  the  near-wake  region  lead  one  to  question  the  use¬ 
fulness  of  an  iterative  procedure  at  present.  It  is  anticipated,  however,  that  once 
improvements  are  made  in  defining  the  displacement  model  over  the  entire  body  length 
and  in  the  wake  region,  an  iterative  procedure  will  be  adopted  again. 

One  further  obstacle  arose  in  defining  the  displacement  body  for  the  3:1  trans- 

2 

verse  cross-section41  model.  Excessive  boundary-layer  growth  in  the  C  K  boundary- 
layer  computation  caused  the  computer  program  to  abort  prematurely.  No  values  for 
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the  effective  displacement  thickness  were  computed  along  the  major  elliptic  axis 
meridian  for  locations  greater  than  81  percent  of  the  body  length.  A  careful  hand¬ 
fairing  was  used  to  define  the  effective  displacement  thickness  along  the  major 
axis  meridian. 


COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL  RESULTS 
All  data  are  presented  in  the  coordinate  system  used  to  experimentally  measure 
the  boundary-layer  flow.  The  coordinate  system,  denoted  x  -  n  -  0,  is  given  in 
Figures  1  and  4.  The  axial  coordinate  x  is  measured  from  the  nose  of  the  body  and 
passes  through  the  center  of  the  elliptic  profile.  The  coordinates  ng  and  9  are 
defined  along  an  axial  cut  normal  to  the  x-axis,  i.e.,  in  the  y-z  plane.  The 
normal  component  n^  is  measured  from  the  model  surface  and  is  normal  to  the  elliptic 
surface.  The  angular  coordinate  9  is  defined  as  the  angle,  in  degrees,  measured 
from  the  z-axis  to  the  line  joining  the  surface  offset  and  elliptic  center. 

PRESSURE  DISTRIBUTION 

The  steady  pressure  was  measured  along  the  stern  surface  using  pressure  taps. 
These  taps  are  located  at  nine  axial  and  five  radial  positions,  for  a  total  of  45 
measurements.  The  pressure  coefficient  C^  is  computed  from  the  measured  pressures 
by  the  relationship 


C 

P 


'— p  p-p 

ro  __  r  ri 

t-ps  ipu 


(4) 


where  p  =  measured  local  static  pressure 
Pq  =  measured  ambient  pressure 
Pt  =  measured  dynamic  total  pressure 
Ps  =  measured  static  pressure 
p  =  mass  density  of  the  fluid 
Uq  =  free-stream  velocity 

The  measured  values  of  the  pressure  coefficients  are  given  in  Table  2  and  com¬ 
pared  in  Figure  6  with  two  analytically-predicted  distributions  of  pressure 
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coefficient.  The  dashed  curve,  denoted  by  potential  flow  theory,  represents  the 
predictions  of  the  XYZ  potential  flow  method  of  Dawson  and  Dean^  before  using  the 
displacement  body  concept.  The  solid  curve  shows  on  the  displacement  body  after 
one  iteration  of  the  displacement  body  procedure.  The  computed  pressure  coefficient 
is 


where  is  the  computed  potential  flow  velocity  on  the  displacement  body  and  U  is 
the  free-stream  velocity,  100  ft/sec  (30.48  m/s). 

Two  results  are  immediately  apparent  from  the  comparisons  given  in  Figure  6. 
First,  the  theory  was  not  able  to  predict  accurately  the  values  of  the  pressure 
coefficient  for  x/L  >  0.93.  At  these  locations,  the  boundary  layer  is  much  thicker 
than  the  body  cross  section  and  theoretical  displacement  thicknesses  were  not 
available  due  to  premature  abortion  of  the  computer  code  calculation  in  the  separa¬ 
tion  region.  Second,  the  predictions  using  the  displacement  body  concept  agree 
more  closely  with  the  measured  values  than  do  the  data  denoted  as  potential  flow. 
After  one  iteration  of  the  displacement  procedure,  overall  agreement  between 
theoretical  and  measured  values  of  the  pressure  coefficient  is  considered  good  even 
though  the  predicted  values  are  slightly  lower  than  the  measured  values.  No  further 
iterations  of  the  displacement  method  have  been  implemented  at  present.  Further 
refinement  of  the  three-dimensional  wake  and  near  wake  region  by  the  displacement 
body  conception  should  improve  the  accuracy  of  the  theoretical  prediction. 

MEAN  VELOCITY  PROFILES 

Mean  velocity  measurements  were  taken  with  an  ”XM  hot-film  sensor  which  was 
stepped  away  from  the  body  in  the  ng  direction.  Measurements  of  velocity  in  the 
axial  x  and  normal  n£  directions,  and  v^,  respectively,  were  taken  with  the 
probe  elements  alined  vertically.  The  sensor  elements  were  rotated  90  degrees  to 
the  horizontal  position  to  measure  the  mean  velocity  Wq  in  the  0  direction.  An 
on-line  computer  was  used  to  collect  data  at  a  sample  rate  of  1024  data  values  in 
8  sec . 
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The  measured  values  of  the  mean  velocity  components  are  listed  in  Tables  3 
through  9  along  with  other  measured  quantities.  Tables  3  through  9  give  the 
measured  data  along  the  0,  67,  80,  83,  86,  87,  and  90-degree  planes,  respectively, 
for  various  axial  locations  along  the  model.  The  velocity  components  are  non- 
dimensionalized  by  the  free-stream  velocity  Uq.  As  shown  in  the  tables,  the  mean 
axial  velocity  is  the  largest  of  the  three  measured  components.  Measured  mean 
velocity  profiles  in  the  x  and  ng  directions  are  shown  in  Figures  7a  through  7c. 

Each  figure  presents  the  profiles  at  various  angular  positions  for  a  particular 
axial  location.  Figure  7a  shows  that  the  mean  axial  velocity  profiles  vary  only 
slightly  with  angular  position  on  the  model  at  x/L  =  0.719.  Also,  the  boundary 
layer  is  thin,  with  an  overall  thickness  of  less  than  1  in.  Little  variation  in 
the  normal  velocity  component  is  noted.  Examining  the  profiles  further  aft  on  the 
model,  the  boundary  layer  thickens  with  increased  angular  position.  Little  varia¬ 
tion  in  profile  occurs  for  angles  less  than  or  equal  to  80  degrees.  However, 
profiles  between  80  and  90  degrees  become  increasingly  fuller  with  increased  angular 
location.  From  repeated  measurements,  the  accuracies  of  the  experimental  measure¬ 
ments  of  >WUo  and  Vn/UQ  are  estimated  to  be  about  0.5  percent  and  1.0  percent, 
respectively. 

Comparisons  of  the  measured  and  predicted  mean  axial  velocity  profiles  are 

shown  in  Figure  8  at  selected  positions  along  the  model.  The  circular  symbols 

represent  the  "X"  hot-film  measurements  and  the  solid  curves  represent  the  theoreti- 

2  6 

cal  results  of  the  C  K  method  using  the  displacement  body  concept.  Calculations 
2 

using  the  C  K  computer  code  were  made  using  the  initial  velocity  profiles  generated 

within  the  computer  code.  Calculations  were  begun  at  1.5  percent  of  the  body  length 

with  the  transition  located  at  3  percent  of  the  body  length.  Use  of  a  limited, 

discrete  set  of  offsets  to  define  the  model  for  computational  purposes  forced  the 

2 

use  of  this  transition  location.  As  shown  in  Figure  8a,  the  C  K  method,  used  with 

and  without  the  displacement  body,  predicted  the  same  profile  at  x/L  =  0.719  and  0 

degrees.  For  the  basic  body  geometry,  prior  to  using  the  displacement  body  concept, 
2 

the  C  K  method  experienced  excessive  boundary- layer  growth  and  aborted  prematurely, 
giving  no  predictions  for  axial  locations  x/L  >  0.81  and  angles  greater  than  80 
degrees.  The  agreement  between  the  computed  and  measured  mean  axial  velocity  pro¬ 
files  is  good  at  x/L  =  0.719  and  0  degrees  where  the  boundary  layer  is  thin. 
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Agreement  is  also  fairly  good  at  x/L  =  0.954  and  0  degrees.  However,  at  x/L  =  0.954 
and  angular  positions  83,  86,  and  90  degrees,  the  measured  axial  velocity  components 
are  smaller  than  the  predicted  components,  with  flow  reversal  predicted  at  83  and  86 
degrees.  Agreement  inside  the  boundary  layer  is  particularly  poor.  Since  the  eddy 
viscosity  model  plays  an  important  role  in  this  region,  it  is  essential  to  examine 
the  eddy  viscosity  model  used  for  computing  the  thick  three-dimensional  stern 
boundary  layer. 

MEASURED  TURBULENCE  CHARACTERISTICS 

The  turbulence  characteristics  of  the  thick  three-dimensional  boundary  layer 
were  measured  using  an  "X"  hot-film  probe.  An  on-line  computer  was  used  to  collect 
data  at  a  sample  rate  of  1024  data  values  in  8  sec.  The  root-mean-square  values  of 
turbulence  velocity  were  recorded  at  each  probe  position  and  the  eddy  viscosity  and 
mixing  length  values  were  computed  from  the  measured  Reynolds  stresses  and  the 
measured  mean  velocity  profiles. 


MEASURED  REYNOLDS  STRESSES 


.1 


The  distribution  of  the  Reynolds  stresses  -u  v  ,  -u  w„,  u  ,  v  ,  and  w„ 

J  x  n  x  0  x  n  0 

represent  the  turbulence  characteristics  in  the  thick  boundary  layer.  The  mean- 

72  7? 

square  turbulent  velocity  fluctuations  u^  in  the  axial  direction  and  v^  in  tne  n£ 


direction,  and  the  Reynolds  stress  -uxvn  were  measured  with  the  "X"  hot-film  probe 
elements  alined  vertically.  The  probe  elements  were  rotated  90  degrees  to  the  hori- 
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zontal  position  to  measure  both  the  turbulent  fluctuation  w^  in  the  0  direction  and 


the  Reynolds  stress  -u^Wq.  Linear  interpolation  was  used  to  approximate  Wq  and 


-u^Wq  at  the  same  off-body  positions  as  the  data  measured  in  the  vertical  direction. 
All  measured  values  of  the  turbulent  fluctuations  and  the  measured  Reynolds  stresses 
are  given  in  Tables  3  through  9. 

The  nondimensionalized  distributions  of  the  measured  turbulent  fluctuations 


J u'2/U  ,  / v'2/U  ,  and  *^w7^ 


_  /U  and  Reynolds  stress  -100  u  v  /U  at  selected 
x  o '  n  o'  do  xno 

locations  along  the  model,  are  shown  in  Figures  9  through  13.  As  can  be  seen 


in  Tables  4  through  8,  the  Reynolds  stress  -u^Wg  is  typically  one  order  of  magnitude 


less  than  the  Reynolds  stress  -u  v 

x  n 


An  exception  to  this  trend  occurs  for  the 


angular  location  of  80  degrees,  where  measured  values  of  -u^Wg  exceed  the  values  of 


-u  v 
x  n 


This  is  the  region  of  predicted  separation  by  the  C  K  computer  code.  The 


measured  distributions  of  -u^Wg  are  not  depicted  graphically. 

The  results  given  in  Figures  9  through  13  and  in  Tables  4  through  8  indicate 
~2 

that  u^  /U  is  the  largest  component  of  turbulent  velocity  fluctuation  and  that  the 

normal  component  v^  /U  is  the  smallest  component.  In  addition,  the  fluctuations 
are  larger  near  the  body's  surface  and  reduce  to  values  near  zero  as  the  edge  of 
the  boundary  layer  is  approached.  At  the  body's  surface,  the  no-slip  boundary  con¬ 
dition  requires  the  velocity  and  turbulent  fluctuations  to  be  zero,  indicating  that 
a  sharp  gradient  exists  in  the  turbulent  fluctuations  at  the  wall.  This  gradient, 
which  becomes  apparent  in  the  measured  data  as  the  boundary  layer  thickens,  is 

evident  at  all  angular  locations  where  x/L  >  0.914.  Similar  trends  have  been  noted 
1  2 

by  Huang  et  al.  '  for  axisymmetric  bodies. 

The  measured  distributions  of  the  Reynolds  stress  -  100  u'v'/U  are  also  shown 

7  x  n  o 

in  Figures  9  through  13.  The  maximum  value  of  this  component  of  Reynolds  stress 
generally  occurs  near  the  body  wall  showing  little  variation  with  location  along  the 
model.  When  the  boundary  layer  is  thin,  the  spatial  resolution  of  the  "X"  hot-film 
probe  may  not  be  fine  enough  to  measure  precisely  the  Reynolds  stress  distributions 

near  the  wall.  The  maximum  value  of  the  -u'v*  Reynolds  stress  occurs  near  the  wall 
for  all  locations  measured  except  x/L  =  0.914  and  0  =  86  degrees. 

^  ^  2  2  2*  2 

A  turbulence  structure  parameter  a,,  where  a.  =  u  v  /q  and  q  =  u'  +  v^  + 
_  1  1  x  n  ^  n  x  n 

>*2  1  2 

Wg  ,  was  investigated  by  Huang  et  al.  ’  for  axisymmetric  bodies.  Huang's  results 

for  axisymmetric  bodies  showed  that  this  parameter  has  a  value  of  0.16  for  0  ng 

0.6  6^  and  that  the  value  of  a^  decreases  toward  the  edge  of  the  boundary  layer. 

The  parameter  6  ,  used  to  normalize  the  distance  from  the  model  n  ,  is  defined  as 

the  distance  from  the  wall  surface  in  the  n  direction  at  which  the  measured  tur- 

e 

-2 

bulent  fluctuation  u^  /U  reaches  the  value  0.01.  Figures  14a  through  14d  show  the 
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range  of  values  of  the  parameter  a^  for  the  three-dimensional  body.  At  most  axial 

positions  for  the  0-  and  67-degree  locations,  the  value  of  a^  is,  approximately, 

0.16  for  n^/5^  0.8.  The  value  of  a^  reduces  to  0.07  at  the  80-degree  plane,  the 

region  of  separation  predicted  by  the  C^K  computer  code.  For  the  remaining  angular 

positions,  the  value  a.  fluctuates  between  0.04  and  0.16.  A  reduction  in  the  value 

1  12 

of  the  parameter  was  also  found  by  Shiloh  et  al.  near  separation  for  an  airfoil 
type  flow. 

The  free-stream  turbulent  velocity  fluctuations  were  not  removed  from  the 
2 

measured  values  of  q  .  The  reduction  in  the  values  of  a^  near  the  edge  of  the 
boundary  layer  may  be  caused,  in  part,  more  by  the  larger  contribution  of  the  free- 

2  > 

stream  turbulence  to  q  than  to  -u'V'  . 

x  n 

EDDY  VISCOSITY  AND  MIXING  LENGTH 

The  values  of  eddy  viscosity  and  mixing  length  are  not  measured  directly,  but 

1  2 

are  obtained,  as  in  the  axisymmetric  case,  *  from  the  measured  values  of  the 

Reynolds  stress  -u'v^  and  the  mean  velocity  gradient  Su^/Sn^.  The  definitions  used 
to  compute  these  quantities  are 


- - T 

u 

x  n 


v 


(6) 


When  the  values  of  wn/u  are  less  than  0.1  and  *  he  value  of  0  is  90  degrees  for  the 

0  x 

present  measurements.  Equation  (6)  may  be  approximated  by 


— t — r 
U  V 

x  n 


3n 

e 


(7) 
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A  spline  curve  is  used  to  fair  the  experimental  data  before  the  velocity  gradient 
is  obtained  numerically. 

The  nondimens ional  distributions  of  the  eddy  viscosity  e/(U^6  *)  determined 
from  the  data  are  shown  in  Figures  15a  through  15d.  The  parameters  and  6^*  are 
defined  as  the  potential  flow  velocity  at  the  edge  of  the  boundary  layer  and  the 
planar  displacement  thickness,  respectively,  for  the  displacement  body.  The  solid 
curve  shown  in  these  figures  is  the  Cebeci  and  Smith^  thin  boundary  layer  formula, 
given  by 


e 

uTS 


0.0168 


1  +  5.5 


n 

e 

T 

r 


(8) 


All  values  of  eddy  viscosity  for  the  3:1  elliptic  model  are  smaller  than  the  experi¬ 
mentally-derived  values  recommended  by  Cebeci  and  Smith  for  thin  boundary  layers. 

The  experimentally-determined  distributions  of  the  nondimensional  mixing  length, 

H  /<$  ,  are  shown  in  Figures  16a  through  16d.  The  solid  curve  in  these  figures 

P  r  13 

represents  the  thin  boundary-layer  model  of  Bradshaw  et  al.  Agreement  between 

theory  and  measurements  is,  at  best,  fair  for  angular  locations  of  0  and  67  degrees; 

for  angular  locations  greater  than  67  degrees,  the  measured  values  of  mixing  length 

are  much  smaller  than  the  predictions. 

1  2 

For  an  axisymmetric  turbulent  boundary  layer,  Huang  et  al.  ’  proposed  a 
turbulence  model  relating  the  mixing  length  to  the  square  root  of  the  entire  tur¬ 
bulence  annulus  area  between  the  body  surface  and  the  edge  of  the  boundary  layer. 

As  seen  in  Figures  14  through  16,  the  values  of  measured  turbulence  intensity,  eddy 
viscosity,  and  intermittency  across  a  turbulent  boundary  layer  decrease  from  a 
maximum  value  at  60  percent  of  the  boundary-layer  thickness  to  zero  at  the  outside 
edge  of  the  boundary  layer.  The  effective  gross  turbulence  area  relevant  to  the 
mixing  length  parameter  is  [ (a+0. 66a) (b+0. 66^)- (a+e^) (b+e^) ] ;  where  and  are 
the  effective  thicknesses  of  the  separation  bubble  (low  turbulence  mixing)  in  the 
direction  of  the  major  and  minor  axes,  a  and  b,  respectively,  of  the  elliptical 
cross-section,  and  6^  and  6^  are  the  boundary-layer  thicknesses  along  the  a  and  b 
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axes.  A  new  mixing  length  model  is  assumed  to  apply  to  a  thick  three-dimensional 
stern  boundary  layer.  The  schematic  representation  of  effective  turbulence  areas, 
as  determined  by  the  areas  between  the  body  surfaces  and  the  contours  of  0.66^  at 
x/L  =  0.81  and  0.95,  are  shown  in  Figure  17.  The  outside  edges  of  the  effective 


turbulence  areas  are  very  close  to  the  contours  of  J  u'^ /U  =  0.04.  Further  out- 

V  x  o 

side  of  these  edges,  turbulence  intensities  reduce  to  0.01  at  the  edge  of  the 
boundary  layer.  The  mixing  length  parameter  is  assumed  to  be  proportional  to  the 
square-root  of  these  effective  turbulence  areas,  e.g.. 


1  ^  v  (a+0. 66  )  (b+0. 6(5,  )-ab  =  A(x) 

3  D 

where  the  value  of  •:  is  assumed  to  be  small  and  will  be  neglected  and  the  value  of 

a  b 

is  zero  since  no  separation  occurs  there.  The  values  of  e and  may  not  be 
negligible  if  the  separation  region  is  so  large  that  the  effective  turbulence  area 
is  reduced  significantly.  However,  in  the  inner  region,  the  conventional  mixing 
length  in  the  wall  region,  Equation  (2),  is  assumed  to  apply.  The  mixing  length  S, 
is  assumed  to  be  the  same  at  the  intersection  of  the  inner  and  the  outer  region, 

y  =  y  in  Equation  (2).  Figures  18a  through  18c  show  the  normalized  mixing  length 

C  12 
distributions  for  three  axisymmetric  bodies  studied  by  Huang  et  al.  ’  These  figures 

show  that  the  measured  values  for  the  three  axisymmetric  models  agree  reasonably 

well;  each  peaking  at  a  value  of  approximately  0.05.  The  values  of  £/A  at  various 

locations  for  the  present  three-dimensional  model  are  shown  in  Figures  18d  through 

I8j.  With  the  exception  of  the  80-degree  angular  location,  values  of  the  non- 

dimensional  mixing  length  remain  fairly  constant  over  the  stern  with  respect  to  both 

angular  and  axial  positions. 

The  data  in  Figure  18  support  the  use  of  a  revised  mixing  length  formulation. 

The  existing  thin  turbulent  boundary-layer  method  can  be  applied  to  the  axisymmetric 

or  three-dimensional  elliptical  body  at  locations  forward  of  where  the  boundary 

layer  thickness  reaches  20  percent  of  the  major  or  minor  axis  value.  Downstream  of 

this  location,  the  apparent  mixing  length  l  may  be  approximated  by  the  thin  flat 

1 3 

boundary  layer  of  Bradshaw  et  al.  as 
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(9) 


l 

-a —  =  1  for  6  <  0.2a  and  6,  <  0.2b 

X/  3  D  — 


2  /[a(x)+0. 65  (x) ] [b(x)+0. 66  (x) ]-a (x)b(x) 

— —  =  - - -  for  §  >  0.2a 

°  >/fa(xTH)+0.  6Sa  (xra)  ]  [b  (Xth)+0.  66b (x^)  ] -a  (x^b  (xra)  or  *  0>2b 

b 

where  x  is  the  axial  location  downstream  of  the  initial  location  of  the  thick  stern 
boundary  layer  x,^.  The  beginning  of  the  thick  stern  boundary  layer  is  selected  as 
the  axial  location  where  the  local  value  of  <5^  or  6^  grows  to  the  value  of  0.2a  or 
0.2b,  respectively  (whichever  occurs  first).  The  new  formulation  can  be  incorporated 
into  existing  axisymmetric  and  three-dimensional  turbulent  boundary-layer  differ¬ 
ential  methods  and  must  be  evaluated  for  a  variety  of  stern  boundary  layers  before 
its  validity  can  be  fully  established. 

CONCLUSIONS 

The  results  of  recent  experimental  investigations  of  the  thick  stern  boundary 
layer  on  a  three-dimensional  body  having  3:1  elliptic  transverse  cross  sections  are 
presented.  Comprehensive  boundary  layer  measurements,  including  mean  and  turbulence 
velocity  profiles  and  static  pressure  distributions  are  given  in  detail. 

An  initial  attempt  has  been  made  at  extending  to  three  dimensions  the  Lighthill 
and  Preston  displacement  body  concept  used  to  treat  the  viscid-inviscid  stern  flow 
interaction  on  axisymmetric  bodies.  The  results  of  this  initial  investigation  indi¬ 
cate  that  the  use  of  the  displacement  model  method  significantly  improves  theoretical 
predictions  of  the  measured  pressure  coefficients  on  the  body  surface.  However, 
agreement  between  measured  and  predicted  pressure  coefficients  remains  poor  in  the 
thick  stern  boundary- layer  region  over  the  last  7  percent  of  the  body.  Theoretical 
predictions  of  the  measured  mean  axial  velocity  profiles  are  satisfactory  in  the 
thin  boundary- layer  region,  but  are  generally  larger  than  the  measured  values  when 
the  boundary  layer  thickens.  Refinements  in  the  present  displacement  body  modeling 
scheme  to  determine  the  effective  displacement  thickness  accurately  over  the  entire 
model  surface  and  wake  may  improve  the  pressure  distribution  predictions  in  the 
thick  stern  boundary  layer. 
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Measured  values  of  eddy  viscosity  and  mixing  length  in  the  thick  stern  boundary 
layer  were  found  to  be  smaller  than  values  which  have  been  proposed  for  thin 
boundary  layers.  Because  eddy  viscosity  and  mixing  length  models  play  an  important 
role  in  boundary-layer  calculations ,  a  modification  of  the  theoretical  mixing  length 
model  is  proposed  which  may  improve  the  prediction  of  the  boundary  layer. 

Further  work  in  this  area  is  needed.  A  larger  data  base  of  experimental  results 
on  a  variety  of  three-dimensional  geometries  will  aid  in  the  development  of  improved 
theoretical  models  to  predict  the  viscid-inviscid  stern  flow  interaction.  The  pro¬ 
posed  new  mixing  length  formulation  must  be  evaluated  further. 
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Figure  7a  -  Nondimensional  Axial  Lengths,  x/L  =  0.719  and  0.810 
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Figure  7  (Continued) 
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Figure  7b  -  Nondimens ional  Axial  Lengths,  x/L  =  0.854  and  0.894 
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Figure  8  -  Computed  and  Measured  Mean  Axial  Velocity 
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Measured  Distributions  of  Reynolds  Stresses  at  Angular  Location  6  =  67  Degrees 
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Figure  10  -  Measured  Distributions  of  Reynolds  Stresses 
at  Angular  Location  0  =  80  Degrees 
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-  Measured  Distributions  of  Reynolds  Stresses 
at  Angular  Location  0  =  83  Degrees 
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Figure  12  -  Measured  Distributions  of  Reynolds  Stresses 


Measured  Distributions  of  Reynolds  Stresses  at  Angular  Location  8  =  87  Degrees 


Figure  13a  -  Nond imensional  Axial  Lengths,  x/L  =  0.719,  0.810,  0.854,  and  0.894 
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Figure  14b  -  Angular  Locations  0  =  80  and  83  Degrees 
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Figure  14c  -  Angular  Locations,  9  =  86  and  87  Degrees 


Figure  15  -  Measured  Distributions  of  Eddy  Viscosity 


Figure  15a  -  Angular  Locations,  0  =  0  and  67  Degrees 
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Figure  15c  -  Angular  Locations,  0  =  36  and  87  Degrees 
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Figure  18  -  Proposed  Similarity  Concept  for  Mixing  Length 
of  Turbulent  Boundary  Layer 
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Figure  18d  -  Present  Model,  0  Degree  Plane 
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Figure  18  (Continued) 


Figure  1 8g  -  Present  Model,  83  Degree  Plane 


Figure  18h  -  Present  Model,  86  Degree  Plane 
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Figure  18  (Continued) 


Figure  1 8 j  -  Present  Model,  90  Degree  Plane 


TABLE  1  -  MODEL  OFFSETS  (INCHES) 


X 

y 

z 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

-0.60000 

0.00000 

-0.69647 

-0.60000 

0.18803 

-0.69365 

-0.60000 

0.37610 

-0.68510 

-0.60000 

0.62683 

-0.66439 

-0.60000 

0.87756 

-0.63207 

-0.60000 

1.14918 

-0.58167 

-0.60000 

1.42081 

-0.51066 

-0.60000 

1.67154 

-0.41788 

-0.60000 

1.88048 

-0.30359 

-0.60000 

2.01838 

-0.18007 

-0.60000 

2.08942 

0.00000 

-1.20000 

0.00000 

-0.99022 

-1.20000 

0.26736 

-0.98621 

-1.20000 

0.53*72 

-0.97405 

-1.20000 

0.89120 

-0.94461 

-1.20000 

1.24768 

-0.89865 

-1.20000 

1.63387 

-0.82700 

-1.20000 

2.02006 

-0.72604 

-1.20000 

2.37654 

-0.59413 

-1.20000 

2.67361 

-0.43163 

-1.20000 

2.86967 

-0.25601 

-1.20000 

2.97067 

0.00000 

-1.80000 

0.00000 

-1.21776 

-1.80000 

0.32880 

-1.21282 

-1.80000 

0.65759 

-1.19787 

-1.80000 

1.09599 

-1.16167 

-1.80000 

1.53438 

-1.10515 

-1.80000 

2.00931 

-1.01703 

-1.80000 

2.48424 

-0.89288 

-1.80000 

2.92263 

-0.73066 

-1.80000 

3.28796 

-0.53081 

-1.80000 

3.52908 

-0.31484 

-1.80000 

3.65329 

0.00000 

-2.40000 

0.00000 

-1.41032 

-2.40000 

0.38079 

-1.40459 

-2.40000 

0.76157 

-t. 38728 

-2.40000 

1.26928 

-1.34536 

-2.40000 

1.77700 

-1.27990 

-2.40000 

2.32702 

-1.17785 

-2.40000 

2.87705 

-1.03406 

X 

y 

z 

-2.40000 

3.38476 

-0.84619 

-2.40000 

3.80785 

-0.61474 

-2.40000 

4.08710 

-0.36463 

-2.40000 

4.23095 

C.  00000 

-3.00000 

0.00000 

-1.57974 

-3.00000 

0.42653 

-1.57333 

-3.00000 

0.85306 

-1.55394 

-3.00000 

1.42177 

-1.50698 

-3.00000 

1 .99047 

-1.43365 

-3.00000 

2.60657 

-1.31934 

-3.00000 

3.22267 

-1.15828 

-3.00000 

3.79138 

-0.94784 

-3.00000 

4.26530 

-0.68859 

-3.00000 

4.57809 

-0.40843 

-3.00000 

4.73922 

0.00000 

-4.20000 

0.00000 

-t. 87052 

-4.20000 

0.50504 

-1.86293 

-4.20000 

1.01008 

-1.83996 

-4.20000 

1.68347 

-1.78436 

-4.20000 

2.35685 

-1.69754 

-4.20000 

3.08635 

-1.56219 

-4.20000 

3.81585 

-1.37149 

-4.20000 

4.48924 

-1.12231 

-4.20000 

5.05040 

-0.81534 

-4.20000 

5.42076 

-0.48361 

-4.20000 

5.61155 

0.00000 

-5.40000 

0.00000 

-2.11461 

-5.40000 

0.57095 

-2.10603 

-5.40000 

1.14189 

-2.08008 

-5.40000 

1.90315 

-2.01721 

-5.40000 

2.66441 

-1.91906 

-5.40000 

3.48911 

-1.76605 

-5.40000 

4.31381 

-1.55046 

-5.40000 

5.07508 

-1.26877 

-5.40000 

5.70946 

-0.92174 

-5.40000 

6.12815 

-0.54672 

-5.40000 

6.34384 

0.00000 

-6.60000 

0.00000 

-2.32285 

-6.60000 

0.62717 

-2.31343 

-6.60000 

1.25434 

-2.28491 

-6.60000 

2.09057 

-2.21586 

-6.60000 

2.92680 

-2.10805 

-6.60000 

3.83271 

-1.93997 

-6.60000 

4.73862 

-1.70314 

-6.60000 

5.57485 

-1.39371 

-6.60000 

6.27170 

-1.01251 

-6.60000 

6.73263 

-0.60056 

-6.60000 

6.96856 

0.00000 

-7.80000 

0.00000 

-2.50114 

-7.80000 

0.67531 

-2.49099 

-7.80000 

1.35062 

-2.46029 
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X 

y 

z 

X 

y 

z 

-7.80000 

2.25103 

-2.38594 

-14.40000 

9.21719 

0.00000 

-7.80000 

3.15144 

-2.26985 

-16.20000 

0.00000 

-3.13374 

-7.80000 

4.12689 

-2.08887 

-16.20000 

0.84611 

-3.12102 

-7.80000 

5.10234 

-1.83387 

-16.20000 

1.69222 

-3.08255 

-7.80000 

6.00275 

-1.50069 

-16.20000 

2.82036 

-2.98939 

-7.80000 

6.75309 

-1.09022 

-16.20000 

3.94851 

-2.84394 

-7.80000 

7.24832 

-0.64665 

-16.20000 

5.17066 

-2.61719 

-7.80000 

7.50343 

0.00000 

-16.20000 

6.39282 

-2.29769 

-9.00000 

0.00000 

-2.65320 

-16.20000 

7.52097 

-1.88024 

-9.00000 

0.71636 

-2.64243 

-16.20000 

8.46109 

-1.36596 

-9.00000 

1.43273 

-2.60987 

-16.20000 

9.08157 

-0.81020 

-9.00000 

2.38788 

-2.53099 

-16.20000 

9.40121 

0.00000 

-9.00000 

3.34303 

-2.40784 

-18.00000 

0.00000 

-3.16563 

-9.00000 

4.37778 

-2.21586 

-18.00000 

0.85472 

-3.15279 

-9.00000 

5.41253 

-1.94536 

-18.00000 

1.70944 

-3.11393 

-9.00000 

6.36768 

-1.59192 

-18.00000 

2.84907 

-3.01982 

-9.00000 

7.16364 

-1.15650 

-18.00000 

3.98870 

-2.87289 

-9.00000 

7.68898 

-0.68596 

-18.00000 

5.22330 

-2.64383 

-9.00000 

7.95960 

0.00000 

-18.00000 

6.45789 

-2.32108 

-10.80000 

0.00000 

-2.83767 

-IB. 00000 

7.59752 

-1.89938 

-10.80000 

0.76617 

-2.82615 

-18.00000 

8.54721 

-1 .37987 

-10.80000 

1.53234 

-2.79132 

-18.00000 

9.17401 

-0.81845 

-10.80000 

2.55390 

-2,70696 

-18.00000 

9.49690 

0.00000 

-10.80000 

3.57546 

-2.57525 

-20.40000 

0.00000 

-3.17543 

-10.8 0000 

4.68215 

-2.36992 

-20.40000 

0.85737 

-3.16254 

-10.80000 

5.78884 

-2.08061 

-20.40000 

1.71473 

-3.12356 

-10.80000 

6.81040 

-1,70260 

-20.40000 

2.85788 

-3.02916 

-10.80000 

7.66170 

-1.23691 

-20.40000 

4.00104 

-2.88178 

-10.80000 

8.22355 

-0.73366 

-20.40000 

5.23945 

-2.65200 

-10.80000 

8.51300 

0.00000 

-20.40000 

6.47787 

-2.32826 

-13.20000 

0.00000 

-3.01212 

-20.40000 

7.62102 

-1.90526 

-13.20000 

0.81327 

-2.99989 

-20.40000 

8.57365 

-1 .38414 

-13.20000 

1 .62654 

-2.96292 

-20.40000 

9.20239 

-0.82098 

-13.20000 

2.71091 

-2.87338 

-20.40000 

9.52628 

0.00000 

-13.20000 

3.79527 

-2.73357 

-22.80000 

0.00000 

-3.17543 

-13.20000 

4.96999 

-2.51561 

-22.80000 

0.85737 

-3.16254 

-13.20000 

6.14472 

-2.20852 

-22.80000 

1.71473 

-3.12356 

-13.20000 

7.22908 

-1.80727 

-22.80000 

2.85788 

-3.02916 

-13.20000 

8.13272 

-1.31295 

-22.80000 

4.00104 

-2.88178 

-13.20000 

8.72912 

-0.77876 

-22.80000 

5.23945 

-2.65200 

-13.20000 

9.03635 

0.00000 

-22.80000 

6.47787 

-2.32826 

-14.40000 

0.00000 

-3.07240 

-22.80000 

7.62102 

-1 .90526 

-14.40000 

0.82955 

-3.05993 

-22.80000 

8.57365 

-1 .38414 

-14.40000 

1.65910 

-3.02222 

-22.80000 

9.20239 

-0.82098 

-14.40000 

2.76516 

-2.93088 

-22.80000 

9.52628 

0.00000 

-14.40000 

3.87122 

-2.78828 

-27.60000 

0.00000 

-3.17543 

-14.40000 

5.06946 

-2.56596 

-27.60000 

0.85737 

-3.16254 

-14.40000 

6.26769 

-2.25272 

-27.60000 

1.71473 

-3.12356 

-14.40000 

7.37376 

-1.84344 

-27.60000 

2.85788 

-3.02916 

-14.40000 

8.29548 

-1.33923 

-27.60000 

4.00104 

-2.88178 

-14.40000 

8.90381 

-0.79434 

-27.60000 

5.23945 

-2.65200 
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> 

I 

t 

I 

\ 

! 

) 

i 


-27.60000 
-27.60000 
-27.60000 
-27.60000 
-27.60000 
-31.20000 
-3) .20000 
-31 .20000 
-31.20000 
-31 .20000 
-31 .20000 
-31 .20000 
-31 .20000 
-31.20000 
-31 .20000 
-31 .20000 
-36.00000 
-36.00000 
-36.00000 
-36.00000 
-36.00000 
-36.00000 
-36.00000 
-36.00000 
-36.00000 
-36.00000 
-'*6. 00000 
"69.60000 
-39.60000 
-39.60000 
-39.60000 
-39.60000 
-39.60000 
-39.60000 
-39.60000 
-39.60000 
-39.60000 
-39.60000 
-44.40000 
-44.40000 
-44,40000 
-44.40000 
-44.40000 
-44.40000 
-44.40000 
-44.40000 
-44.40000 
-44.40000 
-44.40000 
-50.40000 
-50.40000 


6.47787 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 
1.71473 
2.85788 
4.00104 
5.23945 
6.47787 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 
1.71473 
2.85788 
4.00104 
5.23945 
6.47707 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 
1.71473 
2.85788 
4.00)04 
5.23945 
6.47787 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 
1 .7 1473 
2.85788 
4.00104 
5.23945 
6.47787 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 


-2.32826 

-1.90526 

-1.38414 

-0.82098 

0.00000 

-3.17543 

-3.16254 

*3.12356 

-3.02916 

-2.88178 

-2.65200 

-2.32826 

-1.90526 

-1.38414 

-0.82098 

0.00000 

-3.17543 

-3.16254 

-3.12356 

-3.02916 

-2.88178 

-2.65200 

-2.32826 

-1.90526 

-1.38414 

-0.82098 

0.00000 

-3.17543 

-3.16254 


-50.40000 

-50.40000 

-50.40000 

-50.40000 

-50.40000 

-50.40000 

-50.40000 

-50.40000 

-50.40000 

-56.40000 

-56.40000 

-56.40000 

-56.40000 

-56.40000 

-56.40000 

-56.40000 

-56.40000 

-56.40000 

-56.40000 

-56.40000 

-60.00000 

-60.00000 

-60.00000 

-60.00000 

-60.00000 

-60.00000 

-60.00000 

-60.00000 

-60.00000 

-60.00000 

-60.00000 

-63.60000 

-63.60000 

-63.60000 

-63.60000 

-63.60000 

-63.60000 

-63.60000 

-63.60000 

-63.60000 

-63.60000 

-63.60000 

-68.40000 

-68.40000 

-68.40000 

-68.40000 

-68.40000 

-68.40000 

-68.40000 

-68.40000 

-68.40000 


1 .71473 
2.85788 
4.00104 
5.23945 
6.47787 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 
1 .71473 
2.85788 
4.00104 
5.23945 
6.477B7 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 
1 .71473 
2.85788 
4.00104 
5.23945 
6.47787 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 
1.71473 
2.85788 
4.00104 
5.23945 
6.47787 
7.62102 
8.57365 
9.20239 
9.52628 
0.00000 
0.85737 
1.71473 
2.85788 
4.00104 
5.23945 
6.47787 
7.62102 
B. 57365 


-3.12356 

-3.02916 

-2.88178 

-2.65200 

-2.32826 

-1.90526 

-1.38414 

-0.82098 

0.00000 

-3.17543 

-3.16254 

-3.12356 

-3.02916 

-2.88178 

-2.65200 

-2.32826 

-1.90526 

-1.38414 

-0,82090 

0.00000 

-3.17543 

-3.16254 


-3.12356 
-3.02916 
-2.88178 
-2.65200 
-2.32826 
-1.90526 
-1.38414 
-0.32098 
0.00000 
-3.17543 
-3.16254 
-3.12356 
-3.02916 
-2.88178 
-2.65200 
-2.32826 
-1.90526 
-1.38414 
-0.82098 
0.00000 
-3.17543 
-3.16254 
-3.12356 
-3.02916 
-2.88178 
-2.65200 
-2.32826 
-1.90526 
-1.38414 
-0.82098 
0.00000 
-3.17543 
-3.16254 
-3.12356 
-3.02916 
-2.88178 
-2.65200 
-2.32826 
•1 .90526 
-1 .38414 
0.82098 
0.00000 
3.17543 
3.16254 
3.12356 
3.02916 
2.88(78 
2.65200 
.32826 
.90526 
.38414 
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x  y  z 


-68.40000 

9.20239 

-0.82098 

-68.40000 

9.52628 

0.00000 

-73.20000 

0.00000 

-3.17543 

-73.20000 

0.85737 

-3.16254 

-73.20000 

1.71473 

-3.12356 

-73.20000 

2.85788 

-3.02916 

-73.20000 

4.00104 

-2.88178 

-73.20000 

5.23945 

-2.65200 

-73.20000 

6.47787 

-2.32826 

-7 3.20000 

7.62102 

-1.90526 

-73.20000 

8.57365 

-1.38414 

-73.20000 

9.20239 

-0.82098 

-73.20000 

9.52628 

0.00000 

-76.99344 

0.00000 

-3.17158 

-76.99344 

0.85633 

-3.15871 

-76.99344 

1.71265 

-3.11978 

-76.99344 

2.85442 

-3.02550 

-76.99344 

3.99619 

-2.87829 

-76.99344 

5.23311 

-2.64879 

-76.99344 

6.47003 

-2.32544 

-76.99344 

7.61180 

-1.90295 

-76.99344 

8.56327 

-1.38246 

-76.99344 

9.19124 

-0.81999 

-76.99344 

9.51474 

0,00000 

-8t .56868 

0.00000 

-3.13836 

-81.56868 

0.84736 

-3.12562 

-81.56868 

1.69471 

-3.08710 

-81.56868 

2.82452 

-2.99381 

-81.56868 

3.95433 

-2.84814 

-81.56868 

5.17829 

-2.62105 

-81 .56868 

6.40226 

-2.30108 

-81 .56868 

7.53207 

-1.88302 

-81 .56868 

8.47357 

-1.36798 

-81.56868 

9.09497 

-0.81140 

-81 .56868 

9.41508 

0.00000 

-85.22880 

0.00000 

-3.07684 

-85.22880 

0.83075 

-3.06435 

-85.22880 

1.66149 

-3.02658 

-85.22880 

2.76915 

-2.93512 

-85.22880 

3.87682 

-2.79230 

-85.22880 

5.07678 

-2.56967 

-85.22880 

6.27675 

-2.25597 

-85.22880 

7.38441 

-1.84610 

-85.22880 

8.30746 

-1.34116 

-85.22880 

8.91668 

-0.79549 

-85.22880 

9.23051 

0.00000 

-88.88892 

0.00000 

-2.97905 

-88.88892 

0.80434 

-2.96696 

-88.88892 

1.60869 

-2.93039 

-88.88892 

2.68114 

-2.84183 

-88.88892 

3.75360 

-2.70356 

x  y  z 


-88.88892 

4.91543 

-2748799 

-88.88892 

6.07725 

-2.18427 

-88.88892 

7.14971 

-1.78743 

-88.88892 

8.04343 

-1.29854 

-88.88892 

8.63328 

-0.77021 

-88.88892 

8.93714 

0.00000 

-91.17648 

0.00000 

-2.89825 

-91.17648 

0.78253 

-2.88648 

-91.17648 

1.56505 

-2.85091 

-91.1 764F 

2.60842 

-2.76475 

-91.17646 

3.65179 

-2.63023 

-91.17648 

4.78211 

-2.42051 

-91.17648 

5.91242 

-2.12503 

-91.17648 

6.95579 

-1.73895 

-91.17648 

7.82527 

-1.26332 

-91.17648 

8.39912 

-0.74932 

-91.17648 

8.69474 

0.00000 

-94.83660 

0.00000 

-2.73555 

-94.83660 

0.73860 

-2.72445 

-94.83660 

1.47720 

-2.69087 

-94.83660 

2.46199 

-2.60955 

-94.83660 

3.44679 

-2.48258 

-9 4.83660 

4.51366 

-2.28463 

-94.83660 

5.5B052 

-2.00574 

-94.83660 

6.56532 

-1.64133 

-94.83660 

7.38598 

-1.19240 

-94.83660 

7.92762 

-0.70725 

-94.83660 

8.20665 

0.00000 

-97.12416 

0.00000 

-2.61164 

-97.12416 

0.70514 

-2.60105 

-97.12416 

1.41029 

-2.56899 

-97.12416 

2.35048 

-2.49135 

-97.12416 

3.29067 

-2.37013 

-97.12416 

4.30921 

-2.18115 

-97.12416 

5.32775 

-1.91489 

-97.12416 

6.26795 

-1 .56699 

-97.12416 

7.05144 

-1 .13839 

-97.12416 

7.56854 

-0.67522 

-97.12416 

7.83493 

0.00000 

-100.78428 

0.00000 

-2.37435 

-100.78428 

0.64108 

-2.36472 

-100.78428 

1.28215 

-2.33557 

-<00.78428 

2.13692 

-2.26499 

-100.78428 

2.99168 

-2.15478 

-100.78428 

3.91768 

-1.98298 

-100.78428 

4.84368 

-1.74090 

-100.78428 

5.69645 

-1.42461 

-100.78428 

6.41075 

-1 .03496 

-100.78428 

6.88087 

-0.61387 

-100.78428 

7 . 1 2306 

0.00000 

-103.52940 

0.00000 

-2.16142 
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X 

y 

z 

-103.52940 

0.58358 

-2.15265 

-103.52940 

1.16717 

-2.12612 

-103.52940 

1.94528 

-2.06186 

-103.52940 

2.72339 

-1 .96154 

-103.52940 

3.56634 

-1.80514 

-103.52940 

4.40930 

-1 .58478 

-103.52940 

5.18741 

-1 .29685 

-103.52940 

5.83584 

-0.94214 

-103.52940 

6.26380 

-0.55882 

-103.52940 

6.48426 

0.00000 

-105.81696 

0.00000 

-1.95826 

-105.81696 

0.52873 

-1 .95032 

-105.81696 

1.05746 

-1 .92628 

-105.81696 

1 .76244 

-1.86806 

-105.81696 

2.46741 

-1 .77717 

-105.81696 

3.23113 

-1.63547 

-105.81696 

3.99486 

-1.43582 

-105.81696 

4.69983 

-1.17496 

-105.81696 

5.28731 

-0.85359 

-105.81696 

5.67504 

-0.50629 

-105.81696 

5.87479 

0.00000 

-108.10452 

0.00000 

-1 .72912 

-108.10452 

0.46686 

-1.72211 

-108.10452 

0.93373 

-1.70088 

-108.10452 

1.55621 

-1.64948 

-108.10452 

2.17870 

-1.56922 

-108.10452 

2.85305 

-1.44410 

-108.10452 

3.52741 

-1.26781 

-108.10452 

4.14990 

-1 .03747 

-108.10452 

.  4.66863 

-0.75371 

-108.10452 

5.01100 

-0.44705 

-108.10452 

5.18737 

0.00000 

-110.39220 

0.00000 

-1.47123 

-110.39220 

0.39723 

-1.46526 

-110.39220 

0.79447 

-1.44720 

-110.39220 

1.32411 

-1.40347 

-110.39220 

1.85375 

-1 .33518 

-110.39220 

2.42753 

-1.22872 

-110.39220 

3.00131 

-1 .07873 

-110.39220 

3.53096 

-0.88274 

-110.39220 

3.97233 

-0.64130 

-110.39220 

4.26363 

-0.38038 

-110.39220 

4.41370 

0.00000 

-111.60000 

0.00000 

-1.31694 

-111.60000 

0.35557 

-1.31159 

-111.60000 

0.71115 

-1.29543 

-111.60000 

1.18524 

-1.25628 

-111.60000 

1.65934 

-1.19515 

-111.60000 

2.17294 

-1 .09986 

-111.60000 

2.68655 

-0.96559 

-111.60000 

3.16065 

-0.79016 

X 

y 

z 

-111.60000 

3.55573 

-0.57404 

-111.60000 

3.81648 

-0.34048 

-111.60000 

3.95081 

0.00000 

-114.00000 

0.00000 

-0.97745 

-114.00000 

0.26391 

-0.97349 

-114.00000 

0.52783 

-0.96149 

-114.00000 

0.87971 

-0.93243 

-114.00000 

1.23159 

-0.88706 

-114.00000 

1.61280 

-0.81634 

-114.00000 

1 .99401 

-0.71668 

-114.00000 

2.34589 

-0.58647 

-114.00000 

2.63913 

-0.42606 

-114.00000 

2.83266 

-0.25271 

-114.00000 

2.93236 

0.00000 

-114.84000 

0.00000 

-0.84870 

-114.84000 

0.22915 

-0.84526 

-114.84000 

0.45830 

-0.83484 

-114.84000 

0.76383 

-0.80961 

-114.84000 

1.06937 

-0.77022 

-114.84000 

1.40036 

-0.70881 

-114.84000 

1.73136 

-0.62228 

-114.84000 

2.03689 

-0.50922 

-114.84000 

2.29150 

-0.36994 

-114.84000 

2.45955 

-0.21943 

-114.84000 

2.54611 

0.00000 

-116.04000 

0.00000 

-0.66049 

-116.04000 

0.17833 

-0.65781 

-116.04000 

0.35666 

-0.64970 

-116.04000 

0.59444 

-0.63007 

-116.04000 

0.83222 

-0.59941 

-116.04000 

1.08981 

-0.55162 

-116.04000 

1.34740 

-0.48428 

-116.04000 

1.58517 

-0.39629 

-116.04000 

1.78332 

-0.28790 

-116.04000 

1.91410 

-0.17076 

-116.04000 

1.98147 

0.00000 

-116.88000 
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0.12767 
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0.25534 

-0.46513 

-117.36000 

0.42556 

-0.45107 

TABLE  1  (Continued) 


X 

y 

z 

X 

y 

z 

-117.36000 

0.59579 

-0.42912 

-119.76000 

0.00000 

-0.27944 

-117.36000 

0.78020 

-0.39491 

-119.76000 

0.07545 

-0.27830 

-117.36000 

0.96461 

-0.34670 

-119.76000 

0.15090 

-0.27487 

-117.36000 

1.13484 

-0.28371 
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0.25149 

-0.26657 

-117.36000 

1 .27669 

-0.20611 
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0.35209 
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1.37032 
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0.47666 

-0.11917 
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-0.08657 

-1  IB. 44000 

0.20951 
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-120.24000 
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-0.05135 
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0.34918 

-0.37011 
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0.59583 
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0.03539 

-0.13053 
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0.79148 
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0.07077 
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-0.12502 

-118.44000 
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-0.16912 

-120.48000 

0.16513 

-0.11894 

-118.44000 

1.12436 

-0.10031 

-120.48000 

0.21625 

-0.10946 

-118.44000 

1.16394 

0.00000 

-120.48000 

0.26736 

-0.09609 

-118.92000 

0.00000 

-0.35161 

-120.48000 

0.31454 

-0.07864 

-118.92000 

0.09493 

-0.35018 

-120.48000 

0.35386 

-0.05713 

-118.92000 

0.18987 

-0.34586 

-120.48000 

0.37981 

-0.03388 

-1 18.92000 

0.31645 

-0.33541 

-120.48000 

0.39318 

0.00000 

-118.92000 

0.44302 

-0.31909 

-120.72000 

0.00000 

0.00000 

-118.92000 

0.58015 

-0.29365 

-120.72000 

0.00000 

0.00000 

-118.92000 

0.71728 

-0.25780 

-120.72000 

0.00000 

0.00000 

-118.92000 

0.84386 

-0.21096 

-120.72000 

0.00000 

0.00000 

-118.92000 

0.94934 

-0.15326 

-120.72000 

0.00000 

0.00000 

-118.92000 

1.01896 

-0.09091 

-120.72000 

0.00000 

0.00000 

-118.92000 

1.05482 

0.00000 

-120.72000 

0.00000 

0.00000 

-119.52000 

0.00000 

-0.30195 

-120.72000 

0.00000 

0.00000 

-119.52000 

0.08153 

-0.30073 

-120.72000 

0.00000 

0.00000 

-119.52000 

0.16306 

-0.29702 

-120.72000 

0.00000 

0.00000 

-119.52000 

0.27176 

-0.28805 

-120.72000 

0.00000 

0.00000 

-119.52000 

0.38046 

-0.27403 

-119.52000 

0.49822 

-0.25218 

-119.52000 

0.61599 

-0.22140 

-119.52000 

0.72469 

-0.18117 
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TABLE  2  -  MEASURED  PRESSURE  COEFFICIENTS 


x/L 

0,  Angular  Position  (deg) 

0 

45 

67 

80 

90 

0.719 

-0.0442 

-0.0397 

-0.0385 

-0.0421 

-0.0397 

0.810 

-0.0349 

-0.0385 

-0.0385 

-0.0241 

-0.0193 

0.839 

-0.0277 

-0.0277 

-0.0205 

-0.0144 

-0.0012 

0.854 

-0.0169 

-0.0157 

-0.0085 

-0.0012 

+0.0073 

0.879 

+0.0073 

+0.0061 

+0.0097 

+0.0182 

+0.0278 

0.894 

+0.0230 

+0.0218 

+0.0230 

+0.0339 

+0.0448 

0.914 

+0.0375 

+0.0387 

+0.0448 

+0.0557 

+0.0617 

0.934 

+0.0714 

+0.0714 

+0.0787 

+0.0896 

+0.0932 

0.954 

+0.1150 

+0.1126 

+0.1138 

+0.1211 

+0.1235 
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1  NAVSHIPYD  LBEACH/Lib 

2  NAVSHIPYD  MARE 

1  Lib 
1  Code  250 

1  NAVSHIPYD  PUCET/Lib 


1  NCEL/131 

1  NSWC,  White  Oak/Lib 

1  NSWC,  Dahlgren/Lib 

1  NUSC/Lib 


1  NAVSHIPYD  PEARL/ Code  202.32 

1  NAVSEC,  NORVA/ 6660.03,  Blount 

12  DTIC 

1  AF0SR/NAM 

1  AFF0L/FYS,  J.  Olsen 
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Copies 

2  MARAD 

1  Div  of  Ship  R&D 
1  Lib 

1  NASA/HQ/ Lib 

1  NASA/ Ames  Res  Ctr,  Lib 

2  NASA/Langley  Res  Ctr 

1  Lib 

1  D.  Bushnell 
1  NBS/Lib 

1  NSF/Eng  Lib 

1  LC/Sci  &  Tech 

1  DOT/Lib  TAD-491.1 

2  MMA 

1  National  Maritime  Res  Ctr 
1  Lib 

1  U  of  Bridgeport/E.  Uram 

2  U  of  Cal/Dept  Naval  Arch,  Berkeley 

1  Lib 

1  W.  Webster 

2  U  of  Cal.,  San  Diego 

1  A. T.  Ellis 
1  Scripps  Inst  Lib 

5  CIT 

1  Aero  Lib 
1  T.Y.  Wu 
1  A.J.  Acosta 
1  I.  Sabersky 
1  D.  Coles 

1  City  College,  Wave  Hill/Pierson 

1  Catholic  U  of  Aroer/Civil  & 

Mech  Eng 

1  Colorado  State  U/Eng  Res  Ctr 


Copies 

1  U  of  Connecticut/Scottron 

1  Cornell  U/Shen 

1  Florida  Atlantic  U/Tech  Lib 

2  Harvard  U 

1  G.  Carrier 

1  Gordon  McKay  Lib 

1  U  of  Hawaii/Bretschneider 

1  U  of  Illinois/J.  Robertson 

4  U  of  Iowa 

1  Lib 

1  L.  Landweber 

1  J .  Kennedy 
1  V.C.  Patel 

1  Johns  Hopkins  U/Lib 

1  Kansas  State  U/Nesmith 

1  U  of  Kansas/Civil  Eng  Lib 

1  Lehigh  U/Fritz  Eng  Lab  Lib 

4  MIT 

1  Lib 

1  J.R.  Kerwin 
1  P.  Leehey 
1  J.N.  Newman 

2  U  of  Minn/St.  Anthony  Falls 

1  Lib 
1  R.  Arndt 

1  U  of  Mich/NAME/Lib 
1  U  of  Notre  Dame/Eng  Lib 
1  New  York  U/Courant  Inst/Lib 

3  Penn  State 

1  B.R.  Parkin 
1  R.E.  Henderson 
1  ARL  Lib 
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Copies 

1 

Princeton  U/Mellor 

Copies 

1 

Bethlehem  Steel/New  York/Lib 

1 

U  of  Rhode  Island/F.M.  White 

2 

Boeing  Company/Seattle 

1 

Science  Application,  Inc. 

1  Marine  System 

1  P.  Rubber t 

Annapolis,  MD 

C.  von  Kerczek 

1 

Bolt,  Beranek  &  Newman/Lib 

1 

SIT/Lib 

1 

Exxon,  NY/Design  Div/Tank  Dept 

1 

U  of  Texas/Arl  Lib 

1 

Exxon  Math  &  System,  Inc. 

1 

Utah  State  U/Jeppson 

1 

General  Dynamics,  EB/Boatwr ight 

2 

Southwest  Res  Inst 

1 

Flow  Research 

1  Applied  Mech  Rev 

1  Abramson 

1 

Gibbs  &  Cox/Tech  Info 

3 

Stanford  U 

1 

Grumman  Aerospace  Corp/Lib 

1  Eng  Lib 

1  R.  Street,  Dept  Civil  Eng 

5 

Hydronautics 

1 

1  S.J.  Kline,  Dept  Mech  Eng 

Stanford  Res  Inst/Lib 

1  Lib 

1  E.  Miller 

1  V.  Johnson 

1 

U  of  Virginia/Aero  Eng  Dept 

1  C.C.  Hsu 

1  M.  Tulin 

1 

U  of  Washing ton/ Arl  Tech  Lib 

1 

Lockheed,  Sunnyvale/Waic 

2 

VP  I 

1 

Lockheed,  California/Lib 

1  Dept  Mech  Eng 

1  J.  Sc.hetz,  Dept  Aero  & 

1 

Lockheed,  Georgia/Lib 

2 

Ocean  Eng 

Webb  Inst 

1 

McDonnell  Douglas,  Long  Beach 

1  T.  Cebeci 

1  Lib 

1  Ward 

1 

Newport  News  Shipbuilding/Lib 

1 

Woods  Hole/Ocean  Eng 

1 

Nielsen  Eng  &  Res 

1 

Worchester  Pi/Tech  Lib 

1 

Northrop  Corp/Aircraf t  Div 

1 

SNAME/Tech  Lib 

1 

Rand  Corp 

1 

Bell  Aerospace 

1 

Rockwell  International 

1 

Bethlehem  Steel/Sparrows  Point 

1 

1  B.  Ujihara 

Sperry  Rand/Tech  Lib 
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Copies 

Copies 

Code 

1 

Sun  Shipbuilding/Chief  Naval  Arch 

1 

1561 

1 

Robert  Taggart 

1 

1563 

1 

1564 

1 

TRW  Systems  Group/Lib 

1 

1572 

1 

TRACOR 

1 

1606 

1 

United  Technology/East 

1 

1615 

Hartford,  Conn 

1 

1802.1 

2 

Westinghouse  Electric 

1 

1  Q  AO  O 

1  M. 

S.  Macovsky 

X 

XoUZ  •  Z 

1  Gu  1  ino 

1 

1840 

1 

1843 

CENTER  DISTRIBUTION 

1 

19 

Copies 

Code 

Name 

1 

1940 

1 

1500 

W.B.  Morgan 

1 

1942 

1 

1504 

V. J .  Monacella 

1 

l 

1942 

1942 

1 

1507 

D.S.  Cieslowski 

1 

1942 

1 

1508 

R.S.  Rothblum 

10 

5211.1 

1 

152 

W.C.  Lin 

1 

522.1 

1 

1521 

P.  Pien 

1 

1521 

W.  Day 

1 

522.2 

1 

1522 

G.  Dobay 

1 

1522 

M.  Wilson 

1 

154 

J.  McCarthy 

1 

154 

P.  Granville 

1 

1540.1 

B.  Yim 

1 

1540.2 

R.  Cumming 

30 

1542 

T.T.  Huang 

3 

1542 

N.C.  Groves 

1 

1542 

G.S.  Belt 

1 

1542 

Y.T.  Lee 

1 

1542 

M.S.  Chang 

1 

1544 

T.  Brockett 

1 

1544 

R.  Boswell 

1 

1544 

E.  Caster 

1 

1544 

S.  Jessup 

1 

1544 

K.F.  Lin 

1 

1560 

M.  Martin 

Name 

G.  Cox 
W.E.  Smith 
J .  Feldman 

E.  Zarnick 

T.C.  Tai 
R.J.  Furey 

H.  Lugt 

F.  Frenkiel 
J.  Schot 

H.  Haussling 

M.M.  Sevik 

J.T.  Shen 

B.E.  Bowers 
T.M.  Farabee 
F.E.  Geib 
T.C.  Mathews 

Reports  Distribution 
Unclassified  Lib  (C) 
Unclassified  Lib  (A) 
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